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ABSTRACT 



Solar oscillations frequencies show a distinct change 
with solar activity. The changes in frequencies can 
be used to study the time variation of solar struc- 
ture. We discuss constraints on the changes in so- 
lar structure with time as obtained with helioseismic 
data covering the last six years. The frequency varia- 
tions appear to be dominated by changes in the near- 
surface layers rather than by changes in the structure 
of the deeper layers. 
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1. INTRODUCTION 



The study of the Sun's variability has been of in- 
terest to both geophysicists and astrophysicists. For 
geophysicists the interest lies in studying the effect 
on earth's climate. For astrophysicists, the Sun is a 
fairly typical star with the advantage of proximity 
which allows the detailed study of many phenomena 
which arc important to stellar physics, among them 
being magnetic activity. 

The most important change that takes place in the 
Sun on a time-scale shorter that the human life- 
span, and hence, on a time-span that can directly 
affect us, is the solar magnetic cycle with a periodic- 
ity of 11 years. The total energy output of the Sun 
has been precisely measured for about two decades 
now. These measurements show that solar irradiance 
varies with solar cycle, being highest during maxi- 
mum activity (Willson & Hudson 1988). Short-term 
changes on days to months are directly related to the 
evolution of active regions (Frohlich & Pap 1989). In 
this increasingly technological age the effect of the 
solar magnetic cycle is also being felt as the effect 
of magnetic storms which can overload power-lines. 
Thus one of the aims of studying solar variability has 
been to predict solar activity. There is however, no 
reliable way to predict solar activity since there is as 
yet no accepted model of the solar dynamo. 



Till very recently, the only way to study the solar cy- 
cle was to observe changes that take place at the sur- 
face and make models which reproduce the surface 
phenomena. However, the history of studying solar 
rotation has shown us how misleading just reproduc- 
ing surface phenomena can be (see, e.g., the review 
by Libbrecht & Morrow 1991). Helioseismology has 
given us the means to study what happens inside 
the Sun. Solar oscillations frequencies are known to 
change with time, therefore, we now have the abil- 
ity to see changes taking place within the Sun even 
though we cannot yet predict them. We can now 
follow the changes as a function of time and of both 
local and global magnetic activity. Helioseismic data 
thus give strong constraints on the physical mecha- 
nisms that might be postulated to explain solar vari- 
ability. 

Solar oscillations are thought to be linear and adia- 
batic (although adiabaticity breaks down at the out- 
ermost layers, where the time scale of heat transfer is 
smaller than the period of the oscillations). Different 
modes of oscillation are described by three numbers, 
the radial order n which is usually the number of 
nodes in the radial direction. The angular eigen- 
functions are described by spherical harmonics Yj^, 
where / is the degree and m the azimuthal order, and 
—l<rn<l. The degree / is the number of nodal 
lines on the surface, while m is the number of nodes 
along the equator. Details of the properties of solar 
oscillations may be found in Christensen-Dalsgaard 
(1998). 

The most complete information on the oscillations 
is in principle contained in the individual frequen- 
cies u)„im or frequency splittings uj n i m — lu„i between 
modes in the same multiplet, to n i being the average 
multiplet frequency, which for cases of slow rotation 
(as is the case for the Sun) depends only on the struc- 
ture of the Sun. However, it is more usual to param- 
eterize the rotational splittings in terms of splitting 
coefficients according to the formula 

jmax 

= Vnlm = Vnl + ^2 Cbj(n, I) Vf* , (1) 
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where the basis functions are polynomials related to 
the Clebsch-Gordan coefficients (Ritzwoller & Lavely 
1991, Schou et al. 1994). Thus often instead of 
speaking of frequency splittings, one talks of "a co- 
efficients." Because of the structure of the oscilla- 
tion eigenfunctions, the values of odd-order a coeffi- 
cients arc determined by solar rotation, while those 
of even-order coefficients arc determined by a number 
of reasons, such as, deviations from spherical sym- 
metry, large scale flows, magnetic fields, etc. The 
information on all the angular dependences of solar 
structure and dynamics are in the a coefficients. The 
mean frequency v n i depends only on the spherically 
symmetric structure of the Sun. 

In this paper we discuss what we have learned about 
changes that take place in the solar interior because 
of the solar activity cycle. We discuss pre-SOHO 
and pre-GONG results in § 0, the new data used in 
this work are described in §M inversion techniques 
and results are described in § |], a brief description 
of changes in solar dynamics is given in § ||, and we 
state and discuss our conclusions in § |6[ 

2. PRE-SOHO AND PRE-GONG RESULTS 

That solar oscillations frequencies change with time 
has been known since the last solar cycle. Elsworth et 
al. (1990) studied low-degree modes (I = 0, 1, 2, 3) of 
the Sun as obtained by the Birmingham Solar Oscil- 
lations Network (BiSON). They showed that the fre- 
quencies of low degree modes have varied from 1977 
to 1988 in a manner that is correlated with solar ac- 
tivity as measured by the sunspot number. The fre- 
quency variation in that period had a peak-to-peak 
amplitude of 0.46 ± 0.106/iHz. They speculated that 
the change in frequencies could reflect variations in 
the solar radius or in the sound speed in the Sun, and 
these in turn might be due to changes in solar tem- 
perature and/or magnetic fields. The BiSON results 
continue to track solar activity. 

The first results of solar-cycle related changes in in- 
termediate degree modes were reported by Libbrecht 
& Woodard (1990). They showed that measure- 
ments of solar oscillation frequencies in 1986 and 
1988 showed systematic differences of the order of 1 
part in 10,000. The authors averaged the frequency 
differences of modes with degrees of 5 to 50 at given 
frequency intervals, after correcting for differences 
in mode inertia, and showed that the averaged fre- 
quency difference was just a function of frequency, 
and it was a very smooth function of frequency. This 
is a strong indication that the cause of the frequency 
changes is confined to the near-surface layers of the 
Sun. Any perturbation that is localized at some ra- 
dius r m introduces a perturbation to the frequencies 
that is proportional to 

f R dr 

cos(2r m w + <t>), r m = / — (2) 
Jr m c 

(e.ff.. Gouffh 19901. where r m is the acoustical denth 



Table 1. Data sets used 



Set # Start Day 10.7cm Flux' 



GONG sets 
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10-12 
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Jun 
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Sep 
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2000 


167.3 
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Dec 
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2000 


166.1 


MDI sets 
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72.4 
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3. 


1996 


73.2 
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Nov 


28. 


1997 


89.2 
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Apr 
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1998 


108.5 
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Feb 


2. 


1999 


130.7 
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Sep 


7. 


1999 


161.1 


2800 


Sep 


1, 


2000 


175.3 
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Jan 


23, 


2001 


164.4 


3160 


Aug 


27, 


2001 


219.8 



Units of 10~ 22 J s' 1 m- 2 Hz- 1 



of the perturbation, c is the sound speed, u> the angu- 
lar frequency of the oscillation and </> a phase. Hence 
the perturbation contributes to the frequency differ- 
ences a component which is an oscillatory function 
of lo with a 'frequency' / = 2r m which increases with 
the depth of the perturbation. Thus provided they 
have no intrinsic frequency dependence of their own, 
changes in the near surface layers are expected to 
introduce frequency shifts that are slowly varying 
functions of frequency, which is what is seen in the 
Libbrecht & Woodard (1990) data. The authors also 
found changes in the even-order a coefficients and 
speculated that the changes in frequencies and split- 
tings are a result of latitude-dependent changes in 
the structure of the sun which are correlated with 
the Sun's magnetic-activity cycle. 

The Mark-I instrument at the Observatorio del Teide 
have yielded low-degree frequencies for a period of 
fifteen years between 1984 and 1999. These data 
too show that the frequencies track solar activity 
(Jimenez-Reyes et al. 2001). These data also show a 
hysteresis like behavior in the frequency differences 
between the rising and falling phases of Cycle 22 
(Jimenez- Reyes et al. 1998). 

Although oscillation frequencies of low degree modes 
span a long time period, the Big Bear data of Lib- 
brecht, Woodard & Kaufman (1990) and Libbrecht & 
Woodard (1990) are the only widely available inter- 
mediate and high degree data set from the Cycle 22. 
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Figure 1. The frequency differences between MDI 
sets 3160 and 1216 are shown in Panel (a). Panel 
(b) shows the same frequency differences scaled to 
eliminate the effects of mode inertia. Only a few 
error bars are shown for the sake of clarity. 

The situation with Cycle 23 is different, with data 
from both GONG and SOHO, and hence we are now 
in a position to look in detail at changes that take 
place within the Sun as the solar cycle progresses. 



3. NEW DATA 

With the commissioning of the GONG project in 
1995 and the helioseismology instruments on board 
SOHO, such as MDI, GOLF etc., in 1996 we now are 
able to study the time variation of solar oscillations 
frequencies in detail. In this work I shall concentrate 
on results obtained with GONG and MDI. The so- 
lar oscillations frequencies from the GONG project 
were obtained from 108 day time series. The MDI 
data were obtained from 72-day time series. Some 
details of how GONG and MDI projects obtain the 
frequencies are available in Schou et al. (2002). The 
data sets used in this work are listed in Table h The 
heading "Start Day" is the beginning of the 108 day 
observational period for the GONG sets and the 72 
day period for the MDI set. Also listed in the table 
is the average 10.7 cm radio frequency flux over the 
observing period. The 10.7 cm radio frequency flux 
is known to be a good solar activity index. It should 
be noted that only a subset of the available GONG 
and MDI data have been used in this work. 

Figure [l] shows the frequency differences between 
MDI data from late 2001 and mid 1996. Panel (b) of 
the figure shows the frequency differences corrected 



for effect of mode inertia. Modes with low inertia 
are easier to perturb than modes with high inertia. 
The quantity Q n i is the ratio of the mode inertia of a 
mode of degree I order n to that of a mode of degree 
with the same frequency as the mode of degree / or- 
der n (Christensen-Dalsgaard & Berthomieu 1991). 
Scaling by Q n i eliminates the mode-mass (hence de- 
gree) dependence of the frequency differences. We 
see that the scaled frequency differences seem to be 
a function only of frequency. To bring out the trends 
more clearly, we average the scaled frequency dif- 
ferences in groups of 25. These averaged frequency 
differences are shown in in Fig. |^ where we show dif- 
ferent GONG and MDI sets. In the GONG figure 
we have subtracted out the frequencies of set 10-12 
from the other sets. Set 10-12 has the lowest solar 
activity level among the GONG sets used. In the 
MDI figure, we subtract out the frequencies of set 
1216. Again this set has the lowest activity level. 

From Fig. | we see that both GONG and MDI data 
sets behave in similar ways. We also see that the fre- 
quency differences for the different sets seem to be a 
fairly smooth function of frequency. This points to 
the cause of the differences being fairly close to the 
surface of the Sun. We invert the frequency differ- 
ences to check for changes in solar structure. 

The change in frequency tracks solar activity very 
closely. In Fig. |^ we show a plot of the average change 
in scaled frequency for all MDI modes between 2 mHz 
and 3 mHz. The change in frequency was calculated 
with respect to set 1216. Also plotted in the figure 
is the 10.7cm radio flux, and we can see that the 
frequency changes and the radio flux, which is an 
index of solar activity, track each other very closely. 
This behaviour is the same as that found for low 
decree modes by the BiSON and Tenerife groups (see 



4. INVERSIONS FOR SOLAR STRUCTURE 



4.1. Inversion technique 

An inversion for solar structure (e.g. Dziembowski, 
Pamyatnykh & Sienkiewicz 1990; Dappen et al. 1991; 
Antia & Basu 1994; Dziembowski et al. 1994) gen- 
erally proceeds through a linearization of the equa- 
tions of stellar oscillations around a known reference 
model using the variational principle. The differ- 
ences between the structure of the Sun and the ref- 
erence model are then related to the differences in 
the frequencies of the Sun and the model by kernels. 
Non-adiabatic effects and other errors in modeling 
the surface layers give rise to frequency shifts (Cox 
& Kidman 1984; Balmforth 1992) which are not ac- 
counted for by the variational principle. In the ab- 
sence of any reliable formulation, these effects have 
been taken into account in an ad hoc manner by in- 
cluding an arbitrary function of frequency in the vari- 
ational formulation (Dziembowski et al. 1990). Thus 
the fractional change in frequency of a mode can be 
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Figure 2. The averaged, scaled frequency differences between different different GONG sets and GONG set 
10-12 (left panel), and between different MDI sets and MDI set 1216 (right panel). 
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Figure 3. The change in scaled frequencies of dif- 
ferent MDI sets plotted as a function of time (con- 
tinuous line). Also plotted is the 10.7 cm radio flux 
(points). The frequency differences were taken with 
respect to set 1216. 

expressed in terms of fractional changes in the struc- 
ture of the model and a "surface term" . 

When the oscillation equation is linearized — un- 
der the assumption of hydrostatic equilibrium — the 
fractional change in the frequency can be related to 
the fractional changes in two of the functions that 
define the structure of the models. Thus, 

— - l K ^ {r) im dr+ l KpAr) 7tr) dr 



(cf, Dziembowski et al. 1990). Here Su>i is the dif- 
ference in the frequency uii of the ith mode between 
the solar data and a reference model. The function 
c 2 (r) and p(r) are the sound speed and density pro- 
files. The kernels -Fl * 2 p and K l c2 are known func- 
tions of the reference model which relate the changes 
in frequency to the changes in c 2 and p respectively; 
and Ei is the inertia of the mode, normalized by the 
photospheric amplitude of the displacement. The 
term F SUT {, is the "surface term", and results from 
the near-surface errors. The term f sur f also contains 
information of near surface layers that cannot be re- 
solved with the data being inverted. 

Equation (W constitutes the inverse problem that 
must be solved to infer the differences in structure 
between the Sun and the reference model. The in- 
verse problems can be solved by a number of differ- 
ent techniques, such as a Regularized Least Squares 
(RLS) method (see Antia & Basu 1994 for an im- 
plementation of the method), or the method of con- 
structing Subtractive Optimally Localized Averages 
(SOLA) (see Basu et al. 1996; Rabello-Soares et al. 
1999 for implementation). All the results in this work 
were obtained with the SOLA method. The different 
inversion procedures have been tested thoroughly. In 
particular, possible errors introduced in the results 
due to linearization have been studied and found to 
be negligible (Basu et al. 2000). 

For all inversions, we have used Model S of 
Christcnscn-Dalsffaard et al. (1996) as our reference 
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Figure 4- The differences in squared sound speed 
between the Sun and the reference model. 



model. To give a flavour of what the inverted sound- 
speed difference results look like, we show in Fig. ^ 
the sound-speed differences between the Sun and the 
reference model Model S as obtained by inverting 
MDI frequencies obtained from the first 360 days of 
observations. The figure shows that in most parts of 
the solar interior the model is very similar to the real 
Sun, the sound-speed differences being of the order 
of 0.1%. However, there is a noticeable difference 
just below the base of the convection zone (hence 
forth CZ, which is at approximately O.7IR0). The 
solar sound speed is higher in this region. This fea- 
ture is believed to be caused by an excess of helium 
below the base of the CZ in the model. Helium dif- 
fuses below the CZ base because of gravitational set- 
tling. Standard solar models do not have any mixing 
mechanism below the CZ base hence helium builds 
up, increasing the mean molecular weight and hence 
decreasing the sound speed. The base of the CZ in 
the Sun is however a region of large shear, this is 
the location of the tachocline. The shear mixes the 
material below the CZ base and hence the decrease 
in the solar sound speed is less than the decrease in 
the sound speed of the model. 



4.2. Inversion results 

We are mainly interested in seeing how the results 
shown in Fig. ^ change with time. Since the inver- 
sion process is linear, instead of inverting the fre- 
quencies of each set and subtracting the inverted re- 
sults to find the difference, we invert the frequency 
differences. If there are no changes in solar sound 
speed, the inversion results should be zero, within 
limits of error in the inversion results due to errors 
in the frequencies. There seems to be some system- 
atic differences in the high degree modes of the pre- 
and post-recovery MDI sets. Since we are not yet 
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Figure 5. The difference in sound-speed obtained 
by inverting the GONG frequency changes (upper 
panel) and MDI frequency differences (lower panel). 



sure whether the differences are real or are due to 
errors, we truncate the MDI mode set above I = 120 
in order to carry out the inversions. For both GONG 
and MDI data, we restrict the inverted mode set to 
modes with frequencies below 3.5 mHz to avoid large 
surface terms. 

We show the sound speed differences obtained by in- 
verting the frequency differences in Fig. |. These 
inversion results are reliable from about O.IR© to 
O.95R0. The radius range is imposed by the limita- 
tions of the mode sets used in the inversions. We can- 
not get reliable results at smaller radii due to the lack 
of sufficient numbers of low degree modes. We can- 
not get reliable results for radii larger than 0.95Rg, 
due to lack of high degree modes. We see from Fig. || 
that all the inversion results seem to lie within the la- 
error limits. We do not see any systematic change in 
the sound speed differences with time at any radius, 
and hence, conclude that a,t, least within error limits. 
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Figure 6. The position of the base of the solar con- 
vection zone as a function of time and of the 10.7 cm 
radio flux. 



Korzennik (2000) in a preliminary analysis of solar 
cycle related changes in the Sun reported changes 
at the base of the convection zone. Our inversion 
results in Fig. || fail to show any systematic change 
at the CZ base. However, there are better ways to 
look at the CZ base more closely. 

On seeing the results of Eff-Darwich & Korzennik 
(2000), Basu & Antia (2001) investigated the CZ 
base further, in particular, to see if there is any 
change in the position of the base of the CZ, which 
can be determined quite precisely from helioseis- 
mology (see Basu & Antia 1997 for details of the 
method). Basu & Antia (2001) failed to see any 
change in the position of the CZ base. We have ex- 
tended the work to newer data sets here. The posi- 
tion of the CZ base as a function of time and solar 
activity are shown in Fig. |^, and we cannot detect 
any change in the position of the CZ base either as 
a function of time or as a function of solar activ- 
ity. Basu & Antia (2001) did not find any changes 
in the extent of overshoot below the CZ base either. 
The lack of change at the base of the CZ is not alto- 
gether surprising. If we assume that the changes are 
caused by magnetic fields, very high magnetic fields 
to change the structure by adding to the gas pressure 
of that layer. 



there is no temporal change in the structure of the 
Sun. 



4.4. The latitudinal distribution of sound speed 



There is some deviation from the la limit in the re- 
sults of the outer convection zone, however, we can- 
not conclude that the changes are due to changes in 
solar activity. There are two main reasons that can 
cause spurious deviations from la close to the sur- 
face. The first reason is that we have not included 
high degree modes in the inversions. The MDI sets, 
as has been discussed earlier have been restricted to 
modes with I < 120, the GONG sets are limited to 
I < 150. Unfortunately, it is very difficult to de- 
termine the frequencies of very high degree modes 
through global analyses due to the systematic errors 
involved (e.g., Rhodes et al. 1998). Since inversions 
at any given radius is most sensitive to modes with 
their lower turning points at that radius, we are un- 
able to invert the outer regions of the Sun reliably. 
The second reason is correlation of errors. Even if 
errors in the input frequencies are uncorrelated, the 
inversion process correlated the errors in the inver- 
sion results at one radius with errors in inversion 
errors at others (see Howe & Thompson 1996). This 
problem can be severe at large radii (see Rabello- 
Soares et al. 1999). This causes the inversion results 
at large radii to be pulled on one side of the true so- 
lution instead of being randomly distributed on both 
sides of the true solution. The effect is of the order 
of la which is seen in the inversion results in Fig. |5|. 



4.3. The base of the convection zone 
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Figure 7. The average a-i coefficients of different data 
sets plotted as a function of the corresponding Leg- 
endre component of the solar surface magnetic field. 
Both GONG (triangles) and MDI(circles) data are 
shown. 



The seat of the solar dynamo is believed to be close 
to the ba,se of the convection zone. Eff-Darwich & 



Howe, Komm and Hill (1999) have found linear re- 
lations between the even-order a coefficients and the 
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Figure 8. Grey scale map showing the reconstruc- 
tion of the latitudinal dependence of the surface term 
from MDI, multiplied by 100. Overlaid contours 
show the Kitt Peak unsigned magnetic flux with the 
Bq term subtracted; contour spacing is 10G 

Legendre decomposition components of the surface 
magnetic flux up to as, and showed that these re- 
lations could be extended backwards in time to the 
BBSO data from the previous cycle. This can be see 
in Fig. |t] for the coefficient. The coefficients have 
been averaged as per Eq. (5) of Antia et al. (2001), 
which takes into account the effect of mode inertia. 
Antia et al. (2001) that a similar relation holds true 
for all even-order coefficients upto a 14. Coefficients 
of orders higher than 14 were not studies. 

Inversions for the latitudinal distribution of sound 
speed show no changes with time (Antia et al. 2001). 
This leads to the question again whether all the 
changes in the even a coefficients are confined to 
changes in the surface term. This indeed seems to 
be the case. A plot of the surface term obtained 
from the inversions is shown in Fig. ^. Superposed 
on the plot is the surface magnetic field and we can 
see a clear correlation between the two. Thus as in 
the case of the frequencies, we are forced to con- 
clude that changes in the even-order a coefficients 
are caused by changes in the near-surface regions of 
the Sun. 



5. CHANGES IN SOLAR DYNAMICS 



Although all changes in solar structure seems to be 
confined to the outer layers which we are as yet un- 
able to resolve in inversions, changes in solar dynam- 
ics can be detected fairly easily. 

Solar zonal flows show a change in with time (see 
Schou 1999, Howe et al. 2000a, Antia & Basu 2000, 
Antia & Basu 2001, Basu & Antia, this volume). The 
different low- and intermediate-latitude bands mi- 
grate towards the equator, while high latitude bands 
migrate towards the poles (Antia & Basu 2001, Basu 
& Antia. this volume"). Antia & Basu (2000) found 



that these flows penetrate to about O.IR0, similar 
results are found by Howe et al. (2000a). 

Results of changes deeper in the Sun are more con- 
troversial. Howe et al. (2000b) claim to see a 1.3 year 
oscillation in the solar rotation rate at the base of 
the CZ, which Antia & Basu (2000), Corbard et al 
(2001) and Basu & Antia (this volume) fail to find. 

Solar meridional flows change too. The maximum 
speed of the flow shows a decrease with increase in 
activity (Basu & Antia, this volume). In addition 
there are more complex changes on short time-scales 
(Haber et al. 2001). 



6. CONCLUSIONS 

Newer more precise solar oscillation frequencies have 
confirmed that these frequencies change with time. 
The frequencies show an increase with solar activ- 
ity and changes are found to be remarkably well 
correlated with solar activity indices. The change 
in frequency is almost purely a function of the fre- 
quency of the mode, once mode inertia effects are 
taken into account. This points to changes confined 
to the near-surface layers of the Sun as being the 
cause of the observed changes in oscillation frequen- 
cies. We have only looked at low and intermediate 
degree modes in this work. High degree modes ob- 
tained from local area analysis (See Bogart et al. this 
volume) show that changes in the frequencies of high 
degree modes in active regions are not merely func- 
tions of frequency, indicating the possibility of struc- 
tural changes in the solar interior, in particular in the 
near-surface regions that are not resolved properly in 
current inversions. Frequency splittings also change 
with activity. In particular, the even order splittings 
are almost linearly correlated with the corresponding 
component of the solar magnetic field. 

Results obtained by inverting the frequencies ob- 
tained at different epochs show no significant changes 
in the structure of the solar interior, at least up to 
a radius of O.95R0. Data sets used in this work do 
not contain enough high degree modes to allow a re- 
liable inversion of shallower layers. All changes are 
consistent as being due to errors in the frequencies. 
We do not see any changes at the base of the solar 
convection zone. 

The latitudinal distribution of solar sound speed does 
not show any changes with time. The "surface term" 
is however, highly correlated with the latitudinal dis- 
tribution of the solar surface magnetic field. This re- 
sult seems to suggest that the frequencies and split- 
ting coefficients are changed by the direct effects of 
magnetic fields. 

Although solar structure does not appear to show 
changes with time and activity, solar flows do. Both 
zonal and meridional flows show changes with time. 
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